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ABSTRACT 
The  theory  is  postulated  that  the  variational  technique  for 
discovering  the  laws  governing  a  physical  problem,  formulated  within 
a  complex  domain,  is  primarily  a  result  of  requiring  analyticity  in 
the  variational  formulation.   In  support  of  the  theory,  both  in  low 
frequency  network  theory  and  in  antenna  theory,  examples  are  shown  in 
which  the  governing  circuit  laws  are  obtainable  from  functions  which 
do  not  possess  the  stationary  property  simply  by  requiring  analyticity 
in  the  functions.  The  resulting  circuitry  is  the  same  as  that  found 
from  the  conventional  variational  formulation.   In  consequence  of  the 
examples,  for  any  number  of  trial  functions  used  in  setting  up  the 
currents  along  any  number  of  antennas,  it  becomes  evident  that  the 
current  modes  obey  the  laws  of  ordinary  circuit  analysis.  The 
generalized  circuitry  is  thereby  extended  to  any  thin  wire  antenna 
configuration. 
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INTRODUCTION 

There  seems  to  be  an  increasing  trend  in  applying  variational 
techniques  to  the  solution  of  problems  formulated  within  a  complex 
domain.  Basically,  a  complex  parameter  is  introduced  into  a  function 
as  a  control  of  the  perturbation  of  some  quantity.  The  function  is 
then  optimized  by  requiring  the  vanishing  of  the  first  term  of  a 
Taylor's  expansion  about  some  value  of  the  parameter.  The  perturbation 
thus  determined  is  customarily  considered  to  be  the  best  solution  for 
the  trial  functions  used,  provided  the  original  function  possesses  the 
stationary  property,  that  is,  provided  the  original  function  also  has 
a  vanishing  first  term  in  its  Taylor's  expansion  about  the  zero  value 
of  the  parameter. 

Rumsey^  states  that  there  is  no  assurance  that  the  best  approxi- 
mation to  the  true  value  of  a  quantity  is  determined  by  the  variational 
technique  if  the  assumed  distribution  about  the  true  distribution  is 
completely  arbitrary.  He2  further  states  that  his  reaction  theorem 
produces  the  best  value  in  the  sense  that  an  approximate  source  looks 
like  a  true  source  to  any  source  which  can  be  used  as  an  observation. 
The  latter  statement  implies  that  all  approximate  sources  are  required 
to  obey  the  same  physical  laws  which  apply  to  the  true  sources.  Thus, 
it  seems  likely  that  the  best  results  obtained  by  the  variational 
principle  are  those  in  which  the  governing  physical  laws  are  discovered 
and/or  imposed. 


1.  V.  H.  Rumsey,  "Reaction  Concept  in  Electromagnetic  Theory", 
Phys.  Rev.,  Vol.  94-,  pp.  H83-H91;  June,  1954,  p.  H85. 

2.  Ibid. 
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The  existence  of  a  Taylor's  expansion  requires  the  function  to 

be  analytic  in  the  parameters.  This  requirement  is  a  very  stringent 

3 
one.  The  real  and  imaginary  components  are  not  independent.   Indeed, 

the  requirement  of  analyticity  leads  to  the  same  solution  of  a  physical 

problem  as  that  obtained  by  the  conventional  stationary  requirement. 

Hereafter,  the  word  optimize  will  be  used  in  the  sense  of  dis- 
covering the  physical  laws  governing  a  given  problem.  Several  examples, 
both  in  low  frequency  circuit  networks  and  in  antenna  networks,  will 
be  solved  to  illustrate  the  equivalence  of  the  aforementioned  two 
methods  of  optimizing  variational  expressions.  In  this  sense  of 
equivalence,  a  variational  expression  need  not  necessarily  possess  the 
stationary  property. 

In  the  examples,  the  variational  functions  which  are  essentially 
quadric  forms  in  the  parameters  may  be  optimized  by  requiring  analy- 
ticity through  the  existence  of  vanishing  first  derivatives,  whereas, 
the  functions  which  are  essentially  bilinear  forms  in  the  parameters 
and  their  complex  conjugates,  may  be  optimized  by  means  of  requiring 
analyticity  through  the  vanishing  of  the  partial  derivatives  with 
respect  to  the  conjugates  of  the  parameters. 

NOTATION 
For  convenience,  the  shorthand  notation  of  tensor  analysis 
will  be  used.  That  is,  the  summation  symbols  will  be  dropped  and 
summations  will  be  indicated  by  means  of  a  combination  of  superscripts 


3.  P.  M.  Morse,  and  H.  Feschbach,  "Methods  of  Theoretical  Physics", 
McGraw  Hill  Book  Co.,  New  York,  N.  Y. ,  1953,  p.  370. 

U.     Gabriel  Kron,  "Tensor  Analysis  of  Networks",  John  Wiley  and  Sons, 
New  York,  N.  Y. ,  1939. 
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and  subscripts.  The  occurrence  of  a  lower  case  index  as  both  a 
superscript  and  as  a  subscript  implies  summation  over  the  range  of 
the  index,  unless  the  index  is  enclosed  within  a  parenthesis.  A 
capital  index  is  used  to  represent  a  single  value  of  the  range. 
Thus, 

Tff7S 

and  iZ-'iyr-  *s  on^v  a  single  term. 

no 

The  Kronecker  delta, 

J  J 

will  be  used.  Also,  the  impulse  function, 

J°°  %(<+)  *A*  -I  J 


will  be  used.  The  latter  should  not  be  confused  with   £l 

which  is  the  first  variation  of  I.  The  conventional  exp  {Jcvt,    ) 

phasor  notation  will  be  implied. 

LUMPED  NETWORKS 
Some  simple  examples  of  lumped  network  theory  will  be  used 
introductorily . 

Consider  the  network  defined  by, 


2^1***,  i,M~, 


•**  •       a) 


Then,  either  the  quadric  form, 

^    O     a        r*'rJ 


"  +      *~t  >  (2) 


-U- 


or  the  bilinear  form, 


*-±*,jZ'T't 


(3) 


may  be  written.   In  the  bilinear  form,    X    *s  tne  complex 
conjugate  of  -X 

The  real  part  of  W  is  the  time  instantaneous  power  less  the 

5 
time  average  power  loss  per  cycle  in  the  network  '  .  The  real  part 

of  y^~      is  the  time  average  power  loss  and  the  imaginary  part  of 

££?      is  the  peak  reactive  power  of  the  circuit.  Each  has  its  on 

sphere  of  usefulness. 

Given  the  quadric  form  (2),  upon  assuming  the  impedance  matrix, 

// 2/.'//  ,  to  be  symmetric  in  the  subscripts,  the  voltage  laws  for 

the  network  may  be  found  by, 


Likewise,  for  the  bilinear  form  (3), 


f\* 


with    //  i£  .  •}}        not  necessarily  being  symmetric  in  the  subscripts, 
unless  the  reciprocity  theorem  is  implied . 

If,  in  equation  (1),    A\-  ST  yf*.)     £.'       the  network 
becomes  a  two  terminal  one.  The  input  impedance  may  be  written, 

Suppose  it  is  desired  to  find  the  error  produced  in  £ 
if  the    7,       are  replaced  by  approximate  values,  X    "^  C  7f  £T 

being  a  parameter  for  controling  the  relative  size  of  the  approxi- 
mations for  the  currents . 


Rumsey,  op.  cit .  ,  p.  ll^L. 


The  impedance  is  written  as, 

***** '  -^Tx^TTnrP-       '  <7) 

Since    *^i,(o)     «•  J^x,        >  the  error  may  be  determined  by 
expanding  ^j^(g)  in  a  Taylor's  series  about   €»^ 

The  existence  of  a  Taylor's  series  for  a  function  of  a 
complex  variable  requires  the  function  to  be  analytic.  One  way  of 
assuring  analyticity  is  to  specify  the  first  derivative,  for  if  the 

ZL 

first  derivative  exists,  so  do  all  order  derivatives   .   Now, 

7(W~  z**<9*L-&j7jpr    -*   <«    J,   (8) 


or 


'■i««ft'S||r-i**«J  • 


Inspection  of  (9)  reveals  that  perhaps  the  simplest  solution 
is  to  require  the  first  derivative  to  vanish.  Further  inspection 
reveals  that  since, 

cm 


the  derivative  will  vanish  provided, 

z  .  Is-  <  s*  , 

<~£J    -*"       Af    <,  (10) 

Consequently,  the  voltage  law  for  the  network  has  been  determined 
by  requiring  analyticity  at  €■  x&-  O     . 

Proceeding,  the  nth  derivative  is,   , 

6.  Morse  and  Feschbach,  op.  cit.,  p.  374. 


from  which 


Because    £  X    z  %  ?J  appears  as  a  second  order 


(12) 


factor  in  (12),  the  quadric  form  (6)  is  said  to  possess  the  stationary- 
property.  This  is  further  augmented  by  the  fact  that  it  can  be  shown 

that  the  vanishing  of  the  first  derivative  of  an  analytic  function 

7 
produces  a  minimax  in  its  absolute  value.   However,  the  bracketed 

expression  in  (12)  is  the  difference  between  the  impedance  that  would 
be  mathematically  seen  at  the  driving  point,  provided  the  true  currents 
vanished  identically,  and  the  true  impedance  at  the  driving  point.  If 
the  y^     are  completely  arbitrary,  this  difference  may  become  suffi- 
ciently large  to  produce  first  order  errors  in  the  impedance. 

Instead  of  from  the  quadric  form  (2),  the  impedance  may  be  found 
from  the  bilinear  form  (3), 


(13) 


/I*/2 

as  before,  write, 


(14) 


7.  Appendix  A. 
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A  necessary  and  sufficient  condition  that  ^  (?)       De  an 

rf 

analytic         function  of  £.  L4)  is  that  all  order  partial  derivatives 
with  respect  to     ^       vanish,  and  that  the  partial  derivative  with 
respect  to    £.      exists. 
Therefore,   set, 

d  e  <15> 

Ibw, 


(16) 


or, 


)  € 


,^fe^'_  i»W . 


(17) 

As  before,  inspection  reveals  that  if  the  voltage  law  (10)  is 
substituted  into  (17), 

jTST"      Z  r*' 

and  the  derivative  vanishes*  as  a  further  check,  all  higher  order 
derivatives  are  expressable  in  terms  of  the  first  derivative  and 
vanish  when  subjected  to  the  voltage  law.  Consequently,  again  the 
requirement  of  analyticity  at  g  a  O    yields  the  voltage  law  for  the 
network . 

Now,  since  jj?  ($)       is  to  be  analytic, 


*  It  is  assumed  that  the  imaginary  part  of  the  function  does  not 
vanish  identicallv. 
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or, 


J. 


^e         x"i±Hnr         '   -      (19) 


Hence,  unless  the  trivial  requirement  that  7^        be  proportional  to 

_X   Is  imposed,  the  first  derivative  is  nonvanishing  at  ^ts  & 
Since  the  nth  derivative,  subjected  to  (10),  is, 

d.^c»)_  n,/-lm_     f  ±  **  <") 

~^pr  *«l-lxtm1.€l<"J   ~£T~  »       (20) 

the  Taylor's  expansion  about  f?fl     becomes, 

U  (21) 

in  which    %X     r  tf"  •£    ,  and  in  which  the  first  term  within  the 

bracket  is  the  mutual  impedance  between  the   .X    and   ^^ 

referred  to  the  driving  point. 

Although  Ai&h    in  (21)  apparently  is  of  first  order  in  5  J  , 

and  hence  does  not  possess  the  stationary  property,  it  may  be  possible 

to  choose  the  7£         such  that  the  mutual  impedance  is  of  the  same 

order  as  the  self  impedance  ^2*^  .  Therefore,  the  error  in  g.     may 

be  of  the  second  order  even  though  the  bilinear  form  is  not  stationary 

in  the  ordinary  sense.  For  arbitrary  7?^      ,  comparison  of  (21)  with 

(12)  reveals  that  the  variational  error  from  the  bilinear  form  actually 

may  be  smaller  than  the  corresponding  error  from  the  quadric  form. 

The  preceding  discussion  of  the  two  algebraic  forms  frequently 
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encountered  in  physical  problems,  that  is,  the  quadric  form  and  the 
bilinear  form,  very  aptly  illustrates  the  postulated  theory  that, 
optimizing  within  a  complex  domain  is  basically  the  principle  of 
requiring  analyticity,  and  that  its  goal  is  to  discover  the  physical 
lavs  governing  the  subject  under  consideration. 

THE  CYLINDRICAL  ANTENNA 
For  a  perfectly  conducting  thin  cylindrical  antenna,  the  driving 
point  impedance  may  be  found  from  , 


J_ 


in  which ? 


/£   =  axial  coordinate  of  antenna, 
*Y   a  linear  surface  coordinate, 

CL  -     radius  of  antenna, 

/t  »  half  length  of  antenna, 
3>f7/ A*n   -  free  space  wave  number, 
■/Yjy   =  current  distribution  function, 
jFfX)     s  complex  conjugate  of  "f  f kt , 
f\2^       ■  real  part  of  -  -  -   • 


8.  J.  G.  Chaney,  "A  Critical  Study  of  the  Circuit  Concept", 
Jour.  Appl.  Phys.,  Vol  22,  1429-1436,  Dec,  1951. 
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For  thin  antennas,  it  can  be  shown  that  using  the  approximate 

Green's  function,  (y(X--}C  )    ,  yields  the  same  results  as  those 

o 
obtained  by  using  the  exact  Green's  function  '.  That  the  imaginary 

part  of  the  current  product  contributes  nothing  to  ^A        ,  may  be 

verified  by  using  the  conjugate  of  the  assumed  current  product, 

interchanging  the  variables,  and  observing  the  resulting  identity  of 

the  two  integrals  produced  by  the  symmetry  in  the  kernel.  Equation 

(22)  may  be  considered  as  an  extension  of  equation  (13),  and  may  be 

classed  as  a  generalized  bilinear  form. 

Suppose  the  current  function  is  perturbed  while  the  driving 

point  voltage,   /"£   ,  is  kept  fixed;  that  is,  write, 

(23) 
By  formally  writing  the  nth  partial  derivative  of  £       with  respect 

0 

to  C    , 


(24) 

it  is  seen  that  conditions  under  which  the  first  partial  derivative 
vanishes  are  sufficient  for  analyticity.  From, 

*5*  -jg^^fafe.-^^*' 


9 .  C .  T .  Tai ,  "  A  new  interpretation  of  the  integral  equation 
formulation  of  cylindrical  antennas",  IRE  Trans.,  FGAP,  Vol.  AP3; 
pp.  125-127,  July,  1955. 
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write*, 

(26) 
Now,  it  might  be  suspected  that  requiring   z?   to  be  analytic 
at  £  jr  C    would  yield  a  physical  law  for  the  antenna.  Hence,  write, 


(27) 


or, 


fCx) 


G-MtfrW*"- 


(28) 


in  which, 

But,  from  (22), 


+m*&j!*M*w- 


(29) 


Therefore,   from  equations   (28),    (29)  and   (30), 
with      SQd     being  the  impulse  function. 


(30) 


(31) 


I*.  K  «  **  T.  ffo,  (32) 

it  becomes  evident  that  the  Kirchoff  voltage  law  has  been  obtained 


*  It  should  be  remembered  that  €-      must  approach  zero  along  with  €   , 
and  that  (26)  only  is  written  for  explanatory  purposes. 
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for  the  antenna,  namely, 

3o   "* 


J/f  '*  (33) 

For  writing  the  Taylor's  expansion,  determine, 


ale"  " '        * 

and  find, 

SXM 


(34) 


WW      (  *     _  *  ) 


(35) 


in  which,  <•  •  .       •  ^       >»»/^l 


and 


35  f*r*±M%lM  $(*-*)**& 


(36) 


If  the  quadric  form  corresponding  to  (22)  is  vised,  and  if 
analyticity  is  required  by  specifying  some  particular  value  of  the 
derivative,  say  the  simplest  value  apparent,  zero,  examination  of  the 
conditions  for  analyticity  about  £•  =r  o    shows  that  precisely  the 
same  law  (33)  is  obtained.  The  corresponding  Taylor's  expansion  is, 


S  It")  _ 

(37) 


with, 


(33) 
and  as  in  lumped  networks,  it  does  not  necessarily  follow  that  an 
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impedance  for  the  antenna  computed  by  assuming  an  a  priori  current 
distribution  along  the  antenna  in  the  quadric  form  is  any  more  accurate 
than  one  computed  by  assigning  an  a  priori  distribution  along  the 
antenna  in  the  bilinear  form. 

Moreover,  if  an  additional  law  governing  the  current  distribution 
were  discovered  and  imposed  upon  a  solution,  its  accuracy  should  be 
improved.  With  this  in  mind,  return  to  equation  (36). 

Suppose  "f(x)    and  €  7)0(.)    were  not  mutually  independent,  but 
were  components  whose  sum  were  the  true  current  distribution.  There 
would  then  be  no  Taylor's  expansion,  that  is  to  say,  the  expansion 
would  vanish  identically  in  £f  . 

If  the    9/ni[.£(rt  %£(*')]  contributed  nothing  to  the 

integral  in  (36),  the  conjugate  of  Xm7  S^"'*/   could  be  taken 
without  altering  the  value  of  the  integral.  This  would  lead  to, 

>x     3I     J  *  tyj  m  f&  m 

Now,  the  value  of    £*J_fiJ?'    used  in  writing  the  Taylor's 

Hz 

expansion  is, 

If  (39)  is  substituted  for  (36)  in  (40),  and  if    -  J  ° • 

is  evaluated  by  means  of  (31),  it  is  found  that  the  derivative  vanishes, 
and  as  a  consequence,  the  Taylor's  expansion  vanishes  identically. 
Physically,  it  has  been  found  that  the  true  current  modes  obey  the 
reciprocity  theorem.  Consequently,  the  accuracy  of  an  approximation 

-H- 


should  be  increased  by  imposing  the  reciprocity  theorem.  This  imposi- 
tion leads  to  A.&-       in  (35)  subjected  to, 

z,^zt,  -  j-^JfJf  fUL  m  iX(P  l    J  (il) 

The  question  might  arise  as  to  whether  or  not  the  requirement  of 
analytic ity  at  a  nonvanishing  value  of  £,  r  £   would  lead  to  useful 

0 

physical  laws.  If  so,  the  resulting  laws  could  be  used  in  computing 
an  approximate  value  of  2 m    when  neither  T  (x)      nor  ^Pi^J    were 
the  true  current  distribution. 
For  simplicity,  define, 


c      32  ft  ft  3(K)  ft*  '•>  G  (*■*')  <L*  <**', 


and, 


with, 


f.  =  ffci   ,     <?>  -   ?<*)  ■ 

Equation  (23)  becomes, 


(42) 


(A3) 

-15- 


and  equation   (26)   becomes, 

(u) 


< t„*  c  *»0<fr  ef.)-  %  ft,,*  ^  eV<U  =  * 


Equation  (£4)  is  independent  of  £      and  its  solution  is, 
Substituting  into  equation  (4.3),  the  impedance  is  found  to  be, 


^  ro * 


a*. 

The  current  distribution  becomes, 


2'/  *■  2«r 3»"  z»'  <«>> 


r^-i-  fc^y  fS * fc-sJ  S  J. 


U7) 


Now,  equations  (4.6)  and  (47)  reveal  that  the  governing  lav 
applicable  to  circuits  driven  in  parallel  has  been  discovered. 
Thus,  write, 

and  equation  (47)  becomes, 

XL $  ri  P -2  2&  **'J  f(o)  t(*»*' j    pCtiJ    ' 

(49) 

However,   for  arbitrary  -f  {$  and     u'(Jl),  there  is  no  assurance 

i 

that  the  Arrf^  "f(X-)    &( Kl  /   will  contribute  nothing  to  ^?/2_  and 

,7    ,  and  it  is  likely  that  2?,„  x£  £        .To  assure  symmetry 
in  the  subscripts,  the  reciprocity  theorem  will  be  imposed.   Hence, 


-16- 


for  the  final  solution,  the  impedances  are, 

2. 


(50) 
(51) 
(52) 


The  imposition  of  the  reciprocity  theorem  somewhat  may  be 
justified  algebraically  be  recalling  that  it  is  the  impedance  function 
which  has  been  optimized,  and  by  noting  that  the  resulting  formula  (4.6) 
is  symmetric  in  the  subscripts.  Then,  as  far  as  the  equivalent  circuit 
equations  (4.8)  are  concerned  £/x    and  £ZI     could  be  interchanged 
without  altering  the  impedance  £  (  €   )    •  However,  ^^*   and  *•£/ 
would  be  interchanged  in  equation  (49),  resulting  in  a  different 
current  distribution.  If  the  two  current  distributions  were  averaged 
over  the  entire  antenna,  a  resulting  averaging  of  the  corresponding 
driving  point  equations  (48)  would  occur.  This  averaging  would  not 
affect  the  denominator  of  J?0  L^») »  but  i"t  would  affect  the  numerator. 
Thus,  unless  both  -f  (X)  and  <p(K)    are  real,  the  imposition  of  the 
reciprocity  theorem  results  in  a  revision  of  the  current  distribution 
along  the  antenna  with  a  corresponding  revision  of  the  driving  point 
impedance.  This  revision  should  increase  the  accuracy  of  the  solution. 

The  same  governing  circuit  laws  can  be  obtained  from  the  quadric 
form  by  requiring  analyticity  by  specifying  the  existence  of  a  vanishing 
first  derivative.  This  is  the  variational  solution  introduced  by 
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10  31 

Storer   and  subsequently  used  by  Tai   .   However,  neither  writer 

pointed  out  the  physical  significance  ol  the  resulting  formulas. 

Indeed,  Storer   mistakenly   implied  that   /a?0^€J  —  2-0   [  was  being 

made  a  minimum. 

Numerical  results,  based  upon  equations  (4.6),  (4-8),  and  (A9), 
have  been  presented  by  the  writer   for  7  c)c)  ~  S-***  /f(*  ~  /£// 
and  <f()c)=c((t-(M)  e^**'*'  .  Impedances  (50),  (51), 

and  (52)  were  intuitively  substituted  for  the  corresponding  impedances 
of  the  variational  metnod  based  upon  the  quadric  form. 

It  should  be  pointed  out  that  if     L?  L  o)   =  O .         g        and 
£m.     could  be  referred  to  any  value  X,   s  ?C 0     for  which  6f(fC0)  ^£  0 
and  the  current  determined  by  ^PCkj  becomes  parasitically  excited. 

ARBITRARY  NUMBER  OF  CURRENTS  ON  ANTENNA 

It  will  be  shown  that  the  circuit  theory,  obtained  for  the  super- 
position of  two  current  modes  along  an  antenna,  can  be  extended  to  the 
superposition  of  an  arbitrary  number  of  current  modes. 

Let  SS.  be  a  constant  internal  impedance  per  unit  length. 

The  driving  point  impedance  may  be  written, 


10.  J.  E.  Storer,  "Variational  Solution  to  the  Problem  of  the  Symmetri- 
cal Antenna",  Craft  Lab.,  Harvard  Univ.  Cambridge,  Mass.;  Tech. 
Rept.  No.  101,  1950. 

11.  C.  T.  Tai,  "A  Variational  Solution  to  the  Problem  of  Cylindrical 
Antennas",  Stanford  Res.  Inst.,  Stanford,  Cal.,  Tech.  Rept.  No.  12, 
1950. 

12.  Storer,  op.  cit.,  p.  7. 

13.  Appendix  A. 

1A.      J.  G.  Chaney,  "A  Simple  Solution  to  the  Problem  of  the  Cylindrical 

Antenna",  U.S.  Naval  Postgraduate  School,  Monterey,  Cal.,  Tech.  Rept. 
No.  15,  Jan.,  1956. 
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"  (53) 


Let, 


(54) 


J* 

with  the  A  being  parameters  replacing  £  .   Upon  dropping  tne  operator 


Re  and  using  a^.  =  fr(0), 

x 


+  H*  J-l'Li    -Z-ZZ X—  <*(**)****.     (55) 


Let, 


ana 


-*£*  r    t/t .  /^a;  a^;    • 


(57) 


Then,  ,  ^  n,     4, 

Analyticity  will  be  assured  provided,  for  each  r,  the  partial 

— p 
aerivatives  with  respect  to  the  parameters,  A  ,  vanish.  Hence, 


or, 


(59) 


J^fi*dR*~^iJ  "  *°j        <«» 
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and  n  equations  have  been  obtained, 


-t 

Factoring  k    , 


c- tr  T  ~<*     ^     1  „«*> 


(61) 
If,  for  each  t,  the  coefficients  of  A  are  required  to  vanish, 
n  linear  homogeneous  equations  are  obtained  in  the  n  parameters,  namely, 


(62) 
It  will  be  shown  that  equations  (62)  are  linearly  dependent  upon 
a  subset  of  (n  -  l)  linearly  independent  homogeneous  equations. 
For  t  s  1, 

5,  U+  **~   a*>    &U  A"-0.  (63) 

For  t  at  1,  multiply  equations  (63)  by  a^  and  equations  (62)  by  ai# 

After  subtracting, 

(64) 

Thus,  for  each  r,  the  subset  (63)  is  obtained. 

Factoring  A  and  dividing  by  a]_  dquations  (63)  may  be  written  as, 

in  which  I    is  defined  as  a,  \  A   .  A  necessary  and  sufficient 
condition  for  the  solution  of  equations  (65)  for  the  I  is  the  vanishing 
of  the  determinant   /  ^/i^,"-  2-24*1  '     ^^e  vanisning  °f  the  determinant 
may  be  verified  by  noting  that  the  first  row  vanishes  identically. 
Therefore,  equations  (65)  reduce  to 
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C^rEjl'-h  &#  -Z,j)2 J=  ot  C;J  =  2,..-, « 

in  which,  due  to  the  hypothesis  of  n  distinct  current  modes^ 

J»   -a*- 5?  .   .  Thus,  equations  (66)  constitute (n  -  1}  linearly 
independent  equations  sufficient  to  solve  for  the  ratios  I yX,1 
Now,  write  the  n  equations, 


(66) 


(67) 


(68) 


Subtraction  of  the  first  equation  from  each  of  the  other  £i  -  1/ 
equations  results  in  the  system, 

Since  equations  (66)  must  yield  the  ratios  -*•  /«£     for  any 
nonvanishing  I  ,  and  since  the  additional  linearly  independent  equation 
in  (68)  simply  serves  to  assign  a  value  to  I0     in  terms  of  a  parameter 
V  ,  the  ratios  may  be  found  by  solving  the  system  (67). 

Let  l<£-£'  for    /tTzfy--,    <71   ,  Then, 

in  which   ^   is  the  cofactor  of   i?^,*  in  the  matrix:  //iSL   // 
and  in  which  ^  is  the  determinant  \^^^  J  . 

Since  f\  /A   «  ^  /  £»/   ,  the  current  distribution  may  be 
written,  m*        A*W >  f(K) 

XOJ- ^  2.  2~    A     — ^7"T  ' 


(71) 
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Substituting  from  equations  (70)  into  equations  (58),  each  term 
of  (53)  may  be  written, 

**<V  H,   4=,  ^    / 

Substituting     £cn)t+}    for     J//064J    /<3L     5»  £,ld  summin6  OV9r 

r  and  s, 


<?--    V. 


Since,  _,u,    t   _/» Uj 


equation  (72)  becomes, 


/t*-  *    ' 


A 


£  ^ 


rw  v  x  *-  />ir4Vy 


By  equating  from  (74.)  and  (75),  it  is  found  that, 


(72) 


1-1     -**/  (73) 

or  finally, 


\*i   Z£*t  ^  (74) 

From  equations  (71),  if  the  applied  voltage  is   l£  , 


(75) 
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Hence, 


Jfa  r   -*  2    £    A         ,rjw 

^     f\s|    »*/  '^  (76) 


In  consequence  of  equations  (74.)  and  (76),  the  circuit  equations 
(67)  may  be  written  directly  from  lumped  network  theory  and  the  current 
distribution  and  driving  point  impedance  may  be  determined  therefrom. 

The  governing  equations  (67)  can  be  derived  far  more  easily  \>9 
setting  up  the  complex  power  bilinear  form,  perturbing  the  currents, 
using  the  analogue  of  equation  (5),  and  requiring  analyticity  of  the 

The  circuit  equations  also  can  be  found  from  the  analogous  quadric 
form  impedance  function  by  specifying  the  existence  of  a  vanishing  first 
derivative.  The  subsequent  algebra  is  the  same  as  for  the  bilinear  form. 

However,  for  the  bilinear  form,  unless  the  current  modes  are 
natural  modes  and  have  adjusted  themselves  in  accordance  with  the  reci- 
procity theorem,  the  circuit  impedances  in  equations  (67)  will  be 
assymetric  in  the  subscripts.  Since  each  current  mode  should  be  as 
nearly  as  possible,  forced  to  obey  all  the  laws  of  the  true  modes,  the 
reciprocity  theorem  will  be  imposed.  Equations  (57)  become  replaced  by, 


Z    */" 


S^W»« 


(77) 


Perhaps  it  should  be  pointed  out  that,  with  respect  to  the  internal 
impedance,  the  effective  value  of  I  is, 
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(■>)  (78) 

which,  for  r  «  ».  reduces  to, 

7  (79) 

Hence,  the  self  and  mutual  internal  impedances  may  be  computed 
in  the  ordinary  manner  of  coupled  circuits  and  added  to  the  corresponding 
radiation  impedances,  provided  the  bilinear  form  of  the  impedance 
function  is  used. 


COUPLED  ANTENNAS 
The  circuit  laws  governing  an  arbitrary  number  of  coupled  antennas, 
arbitrarily  oriented  with  respect  to  each  other,  and  each  having  an 

arbitrary  number  of  current  modes,  have  been  set  forth  by  the  writer 

15 
in  a  recent  paper   .  The  resulting  circuitry  was  an  intuitive  exten- 
sion of  the  variational  solution  by  Leavis  and  Tai     of  an  arbitrary 
number  of  modes  on  an  arbitrary  number  of  straight  coplanar  antennas, 
each  normal  to  the  line  of  centers.  As  in  the  case  of  previous  varia- 
tional solutions,  the  authors,  Leavis  and  Tai,  did  not  point  out  the 
physical  significance  of  their  solution  with  respect  to  the  generalized 
circuitry. 


15.  J.  G.  Chaney,  "On  coupled  antennas  of  unequal  sizes",  25  May,  1956, 
submitted  for  publication  in  IRE  Trans.  PGAP. 

16.  C.  A.  Leavis  and  C.  T.  Tai,  "A  method  of  analyzing  coupled  antennas 
of  unequal  sizes",  IRE  Trans.  PGAP,  Vol.  AP-4,  pp.  128-132, 
April,  1956. 
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It  is  proposed  to  derive  the  circuit  equations  by  requiring 
analyticity  in  the  parameters  for  the  perturbed  currents  in  the  bilinear 
form  impedance  functions. 

The  driving  point  equations  for  n  antennas  may  be  written, 

£K4t  /"%)*    VK   ,  A,**  /,•"'   %  (80) 


with       , 


*  ^{JxrtM^H^\H 


ind  for  r  ^  s, 


in  which, 


R/  w  x  is  the  distance  between  points, 
X/  \,  X/  x  are  curvilinear  coordinates. 


17.  J.  G.  Chaney,  "A  critical  study  of  the  circuit  concept",  loc.  cit., 
p.  U33. 


Following  Leavis  and  Tai  l   ,    let, 

i'^W*  Aft    (*tohk'  1**'mK>"n*+',-'m»> 

with  e.   r  1,  nL    i  :  n_  +  1.     Thus,   for  each  r,   there  are 
491     _  /y/i    z: /)f  w   M  modes,   and  there  are  a  total  of 

/»    5.    "S^f  /n  -  49      )  modes  altogether. 


(85) 


K*l 


K  =    f  *fo)  . 


Define      Q.       z     f      (O)  .    Then, 


I%)ra*^ 


(86) 


For  ;?_   an  interchange  of  variables  shows  that  it  is  unnecessary  to 
take  the  real  part  of  the  current  product.  For  j? nA   »^e  integration 
can  be  performed  by  dropping  the  real  operator,  repeating  the  integration 
for  the  conjugate  of  the  current  product,  and  averaging  the  two  results. 

Substituting  from  (85)  into  (81)  and  (82),  after  dropping  the 
operator  Re,  ,     -f*i-v      -  M      <*<*) 

(88) 
/A  )   f*J 
Define,         X^    '    °-  ^  <  hj      > 


18.  Leavis  and  Tai,  op.  cit.,  p.  131. 
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t*n   —  ■*,«,«      -  (90) 


and, 


Then,  equations  (87)  and  (88)  may  be  expressed  as, 

Substituting  from  equations  (86),  (89)  and  (91)  into  equations 
(92),  twice  the  complex  power,  #  ,  may  be  obtained  by  summing  over  the 
indices  r  and  s.  Hence, 

2*«W4^B  J^'A\  (93) 


The  governing  circuit  lavs  may  be  determined  by  finding  the  rate 
of  change  of  twice  the  complex  power  with  respect  to  the  conjugate 
currents,   %%    ,  and  requiring  the  antenna  self  and  mutual  impedances 
to  be  analytic  ±n   the  lp    ,  that  is,  by  requiring, 


Since  equation  (93)  is  an  identity, 

2**  *>*«*****  L+    B»      ' 


(95) 
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Substituting  from  equations  (30), 


(96) 


Thus,  there  are   ^V 'W*  -  /      equal  values,  *£  ,  for  the  partial 
derivatives  of  the  right  member  of  (92)  with  respect  to  J.   .  Hence, 
for  each  value  of  p  ,  define  V     *  »  ^  • 

Then,  upon  multiplying  and  dividing  by  £   s   ,   equations  (96) 
become,  h    a  a  k 

(97) 

Equations  (97)  constitute  a  system  for  a  network  of  a    meshes, 

n 

with,  for  each   A.  « 'j  •-•>*' j       ^iC^*,*'  I  meshes  being  driven 

in  parallel.  Therefore,  the  circuit  relationship  has  been  established. 
Let   <-±a  h  be  the  cof actor  of     JL  in  the  matrix 


II 2*  ill  >  and  let       «*  -   /  **  I*/  ' 


(98) 


and, 


(99) 


or, 


"+-*«. 


(ioo) 


Therefore, 

•*  4  &  4s\«  t^t  *  >        (101) 
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and, 


!<*)  =  *    £   / 


1/ 

(102) 


or, 

'  K  (103) 

with, 

/   --2   2     /    -r22    zj  ,     .- 

l-l       J~'   *)+*/*£  I  V*'   4-f      *J  t  (ioa) 


Substituting  from  (103)  into  (80), 
or,  since  (105)  must  hold  for  a  set  of  nonvanishing  1/   , 


(105) 


714.7     ~~    «/t  =='^'  (106) 


and 


4 

Let  J)    be  the  cofactor  of       V  LTj      Jj '  V      // 

J)  rr  J  y+ij  ,     Then, 

2^  -    jf^  ?n  "  J)   *•+  (107) 

Again,  the  accuracy  of  a  solution  should  be  increased  by  requiring 
the  circuit  impedances  to  obey  the  reciprocity  theorem.  That  is,  the 
operator  Re  should  be  applied  to  the  current  products  in  equations  (90) 
and  (91). 

The  dropping  of  the  real  operator  in  equations  (87)  and  (88)  now 

19 
must  be  justified.  For  algebraic  simplicity,  spinor    notation  will 

be  used.  That  is,  a  bar  will  be  placed  over  the  indices  associated 

with  the  conjugate  current  functions  in  equations  (90)  and  (91). 


19.  Gabriel  Kron,  loc,  cit.,  pp.  345-349.' 
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Basically,  this  is  simply  a  means  of  indicating  whether  an  index  is  to 
represent  a  row  or  a  column  in  a  matrix.  Thus,  equations  (97)  may  be 
re-written  as,  A 

in  which  the  dots  indicate  the  missing  columns  in  the  singular  matrices. 

Using  the  conjugate  of  the  currents  in  equations  (87)  and  (88), 
the  complex  power  becomes, 

ZW~  ^Rrf^a*  =   *«*    ******  (1W) 

Differentiating  and  imposing  analyticity, 

*flp.  ***"*>  *«   *   P>«     >  (no) 

and  hence  the  resulting  equations  are, 


<*  rA        —   '  (in) 

in  which  the  dots  indicate  the  missing  rows  in  the  singular  matrices. 

p     G     I, 
Since  the  matrix    //  £  *  Jl      is  the  transpose  of  the  matrix 


\\  Jl.  ~  If         ,  and  since  Jj    T     A  //  multiplies  on  the  1< 


whilst     l  J.  &  0f/         multiplies  on  the  right,  it  follows  that 

equations  (111)  are  identically  equations  (108). 

If  any  assumed  current  should  vanish  at  the  driving  point,  that  is, 

if  some  Cu-  rx—  O  ,  the  corresponding  impedances,  «t  *  )    cou^d  be 

referred  to  any  point  on  the  antenna  where  this  current  is  nonvanishing. 

Since  the  left  member  of  (95)  would  be  independent  of  this  J^_  , 

n 
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the  derivative  of  the  right  member  would  vanish.  Thus,  instead  of 
(96),  for  this  Pr, 

fx,    A.         $4,  (112) 

In  (97),  the  corresponding  r    n.       would  be  zero.  Consequently, 
any  current  which  vanishes  at  the  driving  point  becomes  parasitically 
excited  and  the  generalized  network  theory  is  still  valid  provided  the 
corresponding  K.    is  set  equal  to  zero. 


CONCLUSION 

The  examples  which  have  been  discussed  serve  to  illustrate  the 
theory  that  the  solution  to  a  physical  problem,  set  up  in  a  complex 
domain,  may  be  found  by  requiring  analyticity  in  the  variational  form. 
The  form  need  not  necessarily  possess  the  stationary  property. 

In  consequence,  the  generalized  circuit  theory  based  upon  the 
bilinear  forms  which  are  nonstationary  has  been  extended  to  cover  any 
configuration  of  thin  wire  antennas.  The  solutions  may  possibly  be  more 
accurate  than  solutions  based  upon  quadric  forms  which  are  stationary. 

It  is  believed  that  the  theory  of  optimizing  by  requiring  analyticity 
can  be  carried  over  to  other  fields,  and  it  is  hoped  that  it  will  become 
a  useful  tool  in  solving  physical  problems. 
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APPENDIX  A 


It  wSJJ   be  shown  that  if    W (■&    is  analytic,   setting  £  ^(jj  =r o 
yields  a  minimax  in     /  )jC2)l  % 
Let, 

yyCV  *  DC*,!})  +SI/C*>1). 

A  necessary  condition  that  WCh)   be  analytic  is  that  the 
Cauchy-Riemann  equations  are  satisfied,  that  is, 

Also,  a  necessary  condition  for  an  extremum  in  /  iv(&)  (      is  that, 

«  (2) 

A  necessary  condition  for  equations  (2)  to  hold  simultaneously  is  the 
vanishing  of  the  Jacobian,  that  is, 


=  0 

From  (1)  and  (3), 


(3) 


L  *y   *  (  *yJ       y  u) 


or, 


(5) 
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For  an  analytic  function,  equations  (5)  are  equivalent  to  setting, 

-~T£~a-  (6) 

Write  the  second  derivatives  and  subject  them  to  conditions 


From  (1),  the  LaPlace  equations  follow,  that  is, 


Substituting  from  (11)  into  (3), 


Therefore, 


"** 


(7) 
(8) 


^^Tc  '  (9) 

Substituting  from  (1)  into  (9), 


(10) 


(11) 


(12) 


(13) 


20 
Inequality  (13)  is  the  condition  for  an  extremum  to  be  a  minimax  , 

that  is,  /Ak/  is  neither  a  maximum  nor  a  minimum. 

20.  E.  B.  Wilson,  "Advanced  Calculus",  Ginn  and  Co.,  New  York,  N.  Y., 
pp.  114-13  5,  1912. 
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